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Abstract 



We compute three-point functions of single trace operators in planar A/" = 4 SYM. We 
consider the limit where one of the operators is much smaller than the other two. We find 
a precise match between weak and strong coupling in the Frolov-Tseytlin classical limit for 
a very general class of classical solutions. To achieve this match we clarify the issue of 
back-reaction and identify precisely which three-point functions are captured by a classical 
computation. 
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1 Introduction 

In this paper we consider the computation of planar three-point functions of single trace 
gauge-invariant operators in A/" = 4 SYM. These were studied recently at weak [1, 2, 3, 4] 
and strong coupling [5, 6, 7]. Here, we will take the classical limit of these computations 
and find an exact match. In the study of the two-point functions related to the spectrum 
problem, a similar match was pivotal in establishing a first bridge between weak and strong 
coupling [8, 9, 10, 11, 12]. We hope the same to be true for three-point functions. In this 
introduction we sketch the logical fiow of the paper and present some of our main results. 
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First let us define clearly the object we will be interested in. These are the structure 
constants C123 which appear in the three-point functions 

1 C 
{0,{X,)02{X2)0,{X,)) = ^ |^^^|^,^^,_^3|^^^|^;_^l3_^,|^^^|^3+^,_^, , (1) 

where we normalize the two-point functions as {Oi{x)Oi{0)) = Ixl"^'^^ They are not com- 
pletely unambiguous because we can still multiply Oi by a phase. That would not affect the 
normalization of the two-point function, but it would change Cus by that phase. ^ On the 
other hand the absolute value IC123I is well-defined and is the object we will be computing. 
To remove the dependence on the normalization convention, we will also consider the ratio 



C 



123 



/^ooo 
^123 



(2) 



In this expression C123 is the correlation function between three chiral primary operators 
Of^^ with the same charges as Oi. The three-point function (7^23 is known [13] and are non- 
renormalized. Hence r contains as much information as C123 itself but contains no trivial 
cumbersome combinatoric factors. 

At weak coupling, the leading order planar computation is a purely combinatorial prob- 
lem. At 't Hooft coupling A = there is a large number of operators with the same classical 
dimension. Therefore we need to perform degenerate perturbation theory. That is, we must 
use the eigenvectors of the one-loop dilatation operator that lift the degeneracy [14]. For 
example, for operators made out of two complex scalars Z and X we have 

712 — 1 fields 

^(^)= Yl ^(m, . . . , n^)Tr . . . ZXZ . . . ) (3) 

l<ni<-"<nAr<i^ m-l fields 

where 

V^(ni) = 01 , 

^0(^,^2) = 012 +521021, 
'0(ni, 713, ^3) = 0123 + ^21 0213 + ^32 0132 + 531^23 0312 + ^31^21 0231 + 532^31521 0321 , 

and so on. In these expressions 0^,... = e^^-^i+^^^^2+... g^^d Srs = lot^ 
general we have A! terms. We can think of the scalars X as particles with momenta 
Pr propagating in the Z vacuum. They interact with each other in a purely pairwise way. 
Each scattering event is governed by the two particle S-matrix Srs- These integrable wave 
functions '0(ni, . . . , t^at) go by the name of Bet he wave functions. The momenta Pr are those 
of particles in a finite ring (of length equal to the total number of elementary operators in 
the single trace). Hence they are quantized. The corresponding quantization equations that 
follow from the periodicity of the wave functions are the so-called Bethe equations. One 
often thinks of the states (3) as states in a quantum spin chain [15, 16]. In this picture we 
think of the Z fields as spins up forming a ferromagnetic vacuum while the X fields are spins 
down called magnons. 



^The phase is important for bootstrapping higher-point functions. This goes beyond the scope of the 
current analysis. 
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Figure 1: (a) Three-point function of SU{2) operators at tree level considered in [4]. All contractions 
are such that R-charge is preserved. This is the simplest non-trivial configuration which is not 
extremal. 



The tree-level problem is then trivial to enunciate (although very hard to solve in full 
generality): we need to pick three operators like (3) with different kinds of scalars, as repre- 
sented in figure 1. Then we need to sum over all possible Wick contractions of these three 
operators as represented in the same figure. This problem is complex because the number 
of terms in the Bethe wave functions grows factorially. The final form of the result typically 
involves determinants of the Hessian of the Yang- Yang functional^ and sums over all possible 
ways of partitioning the momenta pj into two sets. They are quite involved; see table 2 of 
[4]. 3 

To make analytic progress at weak coupling it is useful to consider some simplifying 

^The Yang- Yang functional is a functional whose minimization leads to the Bethe equations. 

^Needless to say, the example we just described is not the most general one. A general single trace 
operator can contain any of the six scalars, four covariant derivatives and sixteen fermions of A/" = 4 SYM. 
The operators are then described by a Nested Bethe ansatz and the tools needed to efficiently solve the 
weak coupling combinatorial problem are not yet developed. For example, the main ingredients in the 
combinatorial computation are scalar products between Bethe states. These scalar products are known only 
for SU{2) case. We are currently exploring this problem [17]. 
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Oi = EV'iTr(Z...X...F...) 



Oi ~ U1U2U3 . . . 




03 = ETr(zxy) 




02 = E^2lV {Z...X. ..¥...) 
(a) 



O2 ^ U1U2U3 . . . 



Large J simplification 



Figure 2: The setup considered in this paper for three-point function of SU (3) operators. Two 
operators are large and roughly conjugate to each other while the third operator is small, (a) Exact 
setup. The tree-level contraction of three pure eigenstates [4]. These are very entangled states with 
order (number of X^s andY^s)\ terms. All contractions are such that R-charge is preserved, (h) 
Classical limit. The tree-level contraction of the small operator with the corresponding two large 
coherent states. The coherent states are not entangled at all, so their contraction is trivial. 



limits. In this paper we will consider a classical limit where two of the three operators are 
very large while the third operator is small (compared to the other two). This limit will 
allow us to obtain simple analytic results for the structure constants at weak coupling. Also, 
and this is the main result of the current paper, it will enable the direct comparison with 
the strong coupling results. 

In this paper we consider a generalization of the setup of figure 1 studied in [4] , namely 
we consider SU{?)) operators made out of SU{?)) complex scalars as indicated in (4), see 
figure 2a. 

vacuum excitations notations 
~Oi Z Xandy ={Ji,J2,J3} 

02 Z X and Y F, Z} = { Ji - ji, J2 - J2, J3 + Ja} ^ ' 

03 Z Xand? #{X,F,Z} = Oi,j2,j3} 

The total length of operators Oi and O^, are, respectively, 

J = Ji + J2 + Js and j = ji + j2 + js • 

The operator Os is taken to be small, that is j <C J. For comparison with strong coupling, 
we will focus on the case where Os is a protected chiral primary operator. At weak coupling, 
this condition can be trivially relaxed, see section 3.4. 
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The operators Oi and O2 are taken to be very large and with a lot of excitations such 
that their density is kept fixed in the large length limit. Moreover, we only consider low lying 
excitations around the ferromagnetic vacuum. These are excitations with long wavelength. 
That is, the momenta of the individual excitations of operators Oi and O2 are of order 
1/ J. The energy of the individual excitations is of order 1/ so that the total energy is of 
order 1/J. In short, we will consider the Sutherland classical limit [18], rediscovered in the 
AdS/CFT context in [9]. 

The operator Oi and O2 are taken to be roughly the complex conjugate of each other 
in a sense that will be made precise later. The physical picture is that the large operator 
O2 interacts with the small operator to become the slightly modified final large operator 
Oi. Note that we can not take Oi to be exactly the conjugate of O2. This would lead to a 
vanishing three-point function by R-charge conservation. 

Now, in the classical limit the excitations have long wavelength and can be depicted as 
changing the spin configuration slowly as one goes around the spin chain. Magnon dynamics, 
in this limit, are well approximated by the Landau-Lifshitz field theory [19, 11]. More 
explicitly, we can use coherent states to approximate the exact states (3) or rather their 
Nested Bethe Ansatz generalization. The coherent states mimic the exact Bethe states in 
the following precise sense: when computing with them an average of a classical quantity 
such as the total spin or energy, they yield the same result up to finite size corrections. In 
the spin chain language, 

(^1 1 (^classical 1^1) (^^1 1 ^classical 1 , 1 ^^ /rx 

where ipi denotes the exact Bethe state while (fi denotes the coherent state. More explicitly, 
the coherent states are simply 

Oi{x) = ...(u(") .$)(u("+^) .$) (6) 

where $ = (X^Y^Z) and u^^^ = {u^^\u^^\u^^'^) is a complex vector which is unitary and 
slowly varying, i.e. u^^^ -u^^^ = 1 and u^^^-^^ — u^^^ = 0{1/ J). We can approximate the vector 
u^^^ by a continuous field u(cr) where a = 27rn/J is a continuous parameter taking values 
between and 27r. This continuous field obeys the Landau-Lifshitz equations of motion 
described below. 

Now, the key observation is that a small operator Os is roughly a classical operator 
inserted between operators Oi and (92, see also [2]. Note that this is only true provided 
j ^ J which is the limit we consider. Then in the computation of the Wick contractions 
between the three operators, one can simply replace the large operators by the corresponding 
coherent states. This simplifies the computation enormously. Instead of A^i! terms in the 
Bethe state we have a single term in the coherent state description! The coherent state is 
not entangled at all and the Wick contractions are trivial, see figure 2b. The computation 
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is described in detail in section 3.3 and the result is almost obvious, we find^ 



'-^123 




1 




/^ooo 
'-^123 






Jo 2^ 



(7) 



where Vi = -\/ Ji/J. The notation (7**3 emphasizes that we compute the correlation function 
where operators Oi and O2 are non-protected while Os is protected. As anticipated above, 
the classical limit allowed us to dramatically simplify the exact quantum results. 

We will now switch gears and move to strong coupling. We will show that in a particular 
limit the strong coupling computation precisely reproduces our main result (7). At strong 
coupling each gauge-invariant operator is described by a single string: the large operators are 
described by classical strings, while the small operator is a quantum string. Since we consider 
operators made out of scalars only, the string motion is non-trivial in and is point-like in 
AdS^. The simplest classical string solutions in are point-like strings rotating around one 
of the equators,^ 

U = e'^^v , (8) 

where U = {Xi + iX2, X3 + ^^4, X5 + iXe} and Xi are the sphere embedding coordinates. 
The latter square to 1 so that U is unitary, U • U = 1. In (8), v is a constant unitary vector 
which can be parametrized as 




Ja 



— , J — Ji + J2 + Js ^ 



(9) 



where J = y/X k is the total angular momentum of the point-like string and Ji is the angular 
momentum in the (X2i_i, X2i) plane. The energy of this solution is also equal to = a/A 
In the classical limit we should have k ^ 1. This solution is simply the BMN point-like 
string [20] dual to the BPS operator 

Obps oc Tr [X^'Y^^Z^^) + aU possible permutations . (10) 

A simple generalization of the BMN solution is given by considering solutions of the form 



where k ^ oo while ndrU is held fixed [11, 21]. This is the same [11] as the famous Frolov- 
Tseytlin limit [8] which corresponds to A, J — )• oo with A/J^ <C 1. For example, the energy 
of these solutions E = \/\ n is given by 



E=J 



1 + 



A 

2J2 



'0 



— d^u ■ d^n + 



(12) 



^We are doing quantum mechanics so each state in the Hilbert space is a ray. This translates into a 
U{1) gauge symmetry of individually multiplying the by a phase. Strictly speaking (7) is written in a 

conformal-like gauge introduced in Sec. 3. 3. 

^We will often use the same letter for quantities arising at strong coupling and quantities appearing at 
weak coupling. We do so for quantities that will later be matched between weak and strong coupling. We 
hope this will not raise any confusions. 
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One of the remarkable features of this expansion is that an expansion in powers of A/J^ 
resembles a weak coupling expansion in A. Indeed (12) turns out to coincide precisely with 
the weak coupling spin chain spectrum in the classical limit described by the coherent states 
introduced above! This limit was instrumental in establishing a first bridge between weak 
and strong coupling for non-protected operators in the spectrum problem [8, 9, 10, 11, 12]. 

We will now show that for three-point functions this limit is also extremely useful. It 
allows us to match the strong coupling results with the weak coupling prediction (7) in a 
straightforward way.^ Let us give a fiavor of how it works, the details are presented in 
section 2. At strong coupling, Zarembo computed the three-point function of two large 
classical operators and a small BPS operator [5]. His result leads to 



'-'123 



'-"123 



7^ 



Va/j 



27r 



daUi'Ui^Ui' 

2tT cosh-' KTp 



2^2 



cosh 



where j = ji + j2 + js and 



n 



i + i 



(13) 



(14) 



Here, the U's represent the big classical state (11) and v is given by (9). For a more detailed 
description of the several symbols in (13) see section 2. For large only the first term 
inside the square brackets in (13) contributes, the other two are sub-leading as can be shown 
using the Virasoro constraints. Also, when using (11) we can set r = in u since the 
1/ cosy^^(/^re) factor localizes the integral dit Te — l/tv whereas u changes slowly, see figure 
3. Hence the integral over factorizes and can be done trivially. In this way we obtain 
precisely (7)! This match is valid for any solution in the Frolov-Tseytlin limit. 

In the spectrum problem the generic match of the spectrum problem in the Frolov- 
Tseytlin limit was extremely useful, see [23] for an extensive review. For example, the 
weak/strong coupling match for the spectrum can be derived by showing that both weak 
and strong coupling classical limits are properly described by some algebraic curves [12, 24] 
which become identical in the Frolov-Tseytlin limit. Such algebraic curves were one of the 
main ingredients involved in conjecturing the form of the strong coupling Arutyunov-Frolov- 
Staudacher equations [25] and even the all-loop Beisert-Staudacher equations [26] which 
are important input for the Y-system [27]. Those quantum equations can be interpreted 
as a proper discretization of the classical algebraic curves.'' We hope that the match we 
are observing for structure constants could be equally important in the study of three-point 
functions in A/" = 4 SYM. The next step would be to address the question of the discretization 
of the classical results and then conjecture the all-loop result for the structure constants, at 
least for operators that are asymptotically large. 

In the rest of this paper we will give a more detailed explanation of the arguments above. 
We will then discuss some generalizations and finish with some comments and speculations. 
This paper is organized as indicated in the table of contents. 



^The possibility of obtaining a match for the structure constants by performing this limit was first pointed 
out in [4, 22] 

^This discretization technique lead to a vast number of conjectures for all loop Bethe equations for other 
AdS/CFT systems [32, 28]. 
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Figure 3: The Landau- Lifshitz field theory arises naturally both at weak and strong coupling. How- 
ever, this does not mean that a match of the weak and strong coupling results for the structure 
constants is guaranteed. Indeed, the prescriptions for computing these quantities are a priori very 
different. The weak coupling overlap of the coherent states results in a single integral over the 
position of insertion of the small operator, see (7). The strong coupling computation involves a 
two-dimensional integration of a boundary-to-bulk propagator over the full string worldsheet. In the 
Frolov-Tseytlin limit the integration becomes localized in the slice Tg = and only the integral over 
a survives. In this way the two computations are matched. Note that Tg = Z5 the only point that 
is shared between the Lorentzian and Euclidean solutions; the localization mechanism tell us to take 
a snapshot of the string at precisely this physical point. 

2 Strong coupling 

In this section we give a short review of classical strings in AdS^ x in the Frolov-Tseytlin 
limit. We then present a slight generalization of the three-point function of two large classical 
operators and one small BPS operator considered in [5]. 

In AdS^ we consider trivial point-like solutions sitting in the middle of global AdS with 
global time 

t = I^T. (15) 

On the sphere we consider a generic classical motion. The sphere can be described by three 
complex embedding coordinates 

(f/i,f/2,t/3) = U(a,T) (16) 
constrained by U ■ U = 1. The equations of motion and the Virasoro constraints are given 
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byS [d^da + (a^U • daV) ]U, = and {drU ± d^V) • (a^U ± a^U) = /^^ (see for example 
[29]). An important subclass of classical solutions, which can be mapped to coherent spin 
chain states at weak coupling, can be defined by introducing new coordinates 

U(a, r) = e^-u(a, r) (17) 

and considering the limit [11] 

tv ^ oo with hidrU, d^u held fixed. (18) 
In this limit the equation of motion and the Virasoro constraints reduce to 

2m drU = din + 2u {d^u • d^u) (19) 

and 

u • a^u = . (20) 

The global charges read 

Ja = ^ I — {i^UaUa - lUadrUa) • (21) 

Jo 27r 

Using J = Ji + J2 + -^3, this equation leads to (12) where the energy of the solution is 
E = i<iVX. As explained in detail in the introduction this limit is of great significance for 
establishing a bridge between weak and strong coupling results. 

On the other hand, from a more modern perspective, we now know that the match in the 
Frolov-Tseytlin limit is, to a great extent, a fortunate accident.^ A priori the limit A, J ^ oo 
with X/J^ ^ is not the same as A ^ 0. Indeed, the match of the Frolov-Tseytlin limit 
in the spectrum problem breaks down at some loop order in A/" = 4 SYM.^^ StiU, hopefully 
the first few orders in the computation of structure constants will match as well. Of course, 
given the obvious order of limits issue, a more honest way of interpreting this match is as a 
guiding principle in the pursuit of the all loop result. 



2.1 Holographic three-point functions 

Following Zarembo [5], the holographic three-point function of two heavy operators and a 
light chiral primary operator, see figure 2, is given by 



27r 



^123 — ^ 



da^/XUi'U^'U^' 

271 J COSh-^ tvTp 



2/^2 



cosh HiTp 



(22) 



^We use the signature — h i.e. d^da = —d'^ ^ d^. 

^By taking a certain decoupling limit, it was argued in [30] that the spectrum on both sides of the 
AdS^/CFT/^ correspondence should match at least at the one- loop level. It would be interesting to see if the 
same argument holds for structure constants. 

^^A three loop mismatch was reported in [31]. In more recent examples of AdS/CFT dualities it fails 
already at leading order [32]. However the match can be usually achieved after a universal redefinition of 
the coupling constant. 
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where 



C 



n 



i + i 

2^+2 



(23) 



and V is associated with the point-hke BPS string, see (9), while U are the three complex 
sphere embedding coordinates given in (11). Finally in our limit 



(24) 



is the structure constant for three chiral primaries with the same charges as the classical op- 
erators [13], see appendix B. Note that we are working in Euclidean time = ir and ji, ^2, is 
are respectively the number of X, F, Z fields in operator Oa, such that j = ji + j2 + js- The 
result (22) is a trivial generalization of [5] to a generic SU (3) BPS operator. Moreover, our 
normalization convention is slightly diflFerent. 

Given that we are interested in the Frolov-Tseytlin limit, where ^ oo while K.drU is 
held fixed, the integral over in (22) can be easily computed because only the first term in 
square brackets contributes at leading order, while the remaining two are subleading, as can 
be shown using the Virasoro contraint (20): 



2drrU • drrU 



COSh^ HiTe " COSh^ tvTe " " COSh^ tvTe 

Then, using Ui{a^ r^) = e'^^^Ui{(j^ r^) from (11) in (22), we notice that we can set Tg = in Ui 
because the factor 1/ cosh-^^^(/^re) localizes the integral around Tg ^ 1//^, while u is a slowly 
changing variable. Hence, we can evaluate the integral over as (here T = /^Tg) 



dT 



o(jl+j2-j3)T 



2^- 



cos 



1 



Therefore, the three-point function (22) takes the simple form 

27r 



"-"123 



J- 



(ii+i2)! JiWh^- 



Jl!j2! 




da 
2^ 



ui? ul? 



(25) 



Moreover, if we use the equation that relates IZ and C (23), we find a simple expression for 
the ratio CrisVQ: 



OO . 

123- 



27r 



^123 
/^ooo 
'-^123 



'Ui^ vi^ vi^ 



da 
2^ 



(26) 



Te=0 



3 Weak coupling 

In this section we will describe the weak coupling computation of the structure constants in 
the limit where two of the operators are large classical operators whereas the third is a small 
operator. We start by a very representative example and move to the general treatment later 
on. 
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Figure 4- Two Bethe roots configurations for a folded string with a = 1/3 and J = 42,168; 
represented by the larger and smaller bullets respectively. The horizontal and vertical axes are the 
real and imaginary part of u/J, where u is the rapidity of an excitation, related to its momentum 
by u — ^ cot I . The case J — 42 corresponds to the maximum number of roots that we used to 
compute r with the formulas of [4]- Even though this number of roots was not very large, we see 
that they lie nicely along the cuts formed by the clearly classical configuration J = 168. 



3.1 SU{2) Folded string warm-up 

Before moving to the general setup, we provide in this section a most convincing evidence 
for the vahdity of the weak/strong couphng match. We consider the correlation function 
between the following three operators. For the small operator Os we take the BPS operator 

Os = 2 Tt{ZZXX) + Tt{ZXZX) . (27) 

That is we consider ji = js = 2 and ^2 = 0. For the classical states Oi and O2 we will 
take operators dual to the folded string with unit mode number. The folded string in the 
Frolov-Tseytlin limit is given by [8, 11]-^-^ 

/ X 2(l-g)K(g)2 /2K((?) , \ _ , , 2qK(qf re /2K((?) , \ , 

Ui[a,Te) = e -2 K^snl— — cr|(?l , U3[a,Te) = e . ^n ( — — cr|(? 1 , (28) 
for filling fraction^^ 

„._A^ = I_|(£l ,29) 

J1 + J3 K(g) 

We are in the SU{2) setup of [4] which amounts to setting J2 = 0. We use Mathematica's 
conventions where E(g) = EllipticE[q] , K(g) = EllipticK[q] , sn(x|^) = JacobiSN [x,y] 
and dn{x\y) = JacobiDN [x,y] . 



Coherent states and folded strings have been studied in the past in connection with string splitting and 
joining in [33]. 

-^^For completeness, the energy of this solution is given by ^/Xk = E = J -\-2K{q)(E{q) — (1 — q)K{q)) / ir'^ J . 
We will not make use of this expression. 
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Plugging these expressions into (26) we obtain 
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r = 



7rg(l-g)2Fi(|,|;: 
4a{l-a)K{q) 



2;g) 



(31) 



For different filling fractions a we find q from (29), plug it in this expression and get a 
number. For example, for a = 1/4 we find 



Several such predictions are represented by the gray dashed lines in figure 6. 

The claim of the current paper is that these numbers can be matched with the weak 
coupling structure functions. 

We will first check this claim by comparing (31) with the exact tree-level results of [4]. 
The data required to compute Cus from the results of that paper are the positions of the 
Bethe roots of the three operators. For a folded string the Bet he roots are distributed along 
two symmetric cuts, see figure 4. 

Consider first the operator The total number of Bethe roots is the number of X fields 
which is equal to Ji — 2. In each cut we have (Ji — 2)/2 roots. For Oi we have Ji scalars X. 
Hence, operator Oi is parametrized by Ji Bethe roots, two more roots than O2. We want 
to consider an operator Oi which is roughly the complex conjugate of So, where do we 
put the extra two roots of 

Let us recall what the three natural options are when it comes to adding extra roots to 
a classical configuration (these are depicted in figure 5). 

(A) One option is to put extra roots at infinity. Roots at infinity correspond to acting with 
a global symmetry generator on the classical configuration, see e.g. [34, 26]. In this 
case, putting the two roots at infinity would amount to considering the Oi spin chain 
state to be obtained from the O2 state by acting twice with the lowering operator S~ . 
This is the global symmetry operator which converts a Z field into an X. 

(B) Another option is to add roots at finite values outside the classical cuts. Putting roots 
outside the cuts at finite values can be interpreted as considering quantum fiuctuations 
around the classical solution [35]. 

(C) The last option is to add more roots to the classical cuts which already exist. The 
number of roots on each of these classical cuts corresponds to the action variables of the 
classical solution [36] . Hence adding roots to already open cuts leads to a new classical 
state with slightly larger magnitude for the same excited action variables. 

^^For a general BPS operator Os in the SU{2) sector, that is with j2 = 0, we find 



r ~ 0.856897. 



(32) 



r = 



C-^) qnl^ (1 - 9)^3/2 ^p^ (i^ i±i. 2±n . g) 



(30) 



where ji is an even number so that when O2 is a zero- momentum state, so is 0\. 



13 



(A) (B) (C) 

Figure 5: Three different types of deformations of a classical state: (A) A global symmetry trans- 
formation. (B) A quantum fluctuation. (C) A fluctuation in an already existing cut (i.e. adding a 
'^zero mode'^). 

Naively it might seem like it doesn't matter where we add the extra roots and all of these 
options should correspond to Oi This is not the case. As will be explained later, 

there is a sense in which only the last option corresponds to Oi So, let us consider 

that case and add one root to each of the cuts of (92, which has (Ji — 2)/2 roots per cut. 

We will now provide a numerical check of the agreement between the weak and strong 
coupling results in the Frolov-Tseytlin limit. We will discuss the general analytic match for 
general BPS states and general classical solutions later. For the numerics we need to find 
the solutions to Bethe equations with two symmetric cuts for both Oi and O2. This is very 
simple to do even for a very large number of roots (see for example [37]). The obstacle 
is that we then need to plug these solutions into the expressions of [4] which involve sums 
over partitions of Bethe roots, computations of sub-determinants and so on. In practice this 
means that the doable numbers of roots is not very large. This is not a serious problem: 
what we do is compute the three-point function for as many roots as we can and fit the data 
to estimate the large-length asymptotics. 

For example, for a = 1/4 we find 

data={{8, 1.27416239226}, {16, 0.976868082293}, {24, 0.919415218869}, 
{32, 0.897060720354}, {40, 0.885697954062}, {48, 0.878999431990}, 

{56, 0.874655453369}} 

where the first entry is the length of operator Oi and the second entry is r = | C'^'s /C'm I 
computed from [4], see also appendix A. Performing a fit as 

Fit [data, J/J^Range [7] , J] 

we find the numerical prediction^^ 

Tnum ^ 0.856889 (33) 

^^We can also use the numerical data to predict the 1/J finite size corrections. For example, for filling 
fraction a = 1/4 we find r = 0.856889 + 0.565177/ J + . . . . In principle the first 1/J correction should be 
matched with the (Frolov-Tseytlin expansion of the) one loop computation at strong coupling. It would be 
extremely interesting to try this. 
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10 50 100 500 1000 



Figure 6: Fits of the numerical data (black curves) compared to the analytical prediction (dashed 
gray lines) for different values of the filling fraction of operator Oi: a — 1/3, 1/4, 1/5, 1/6; from 
bottom to top. 



a 


Prediction 


Extrapolation 


Error 


1 

f 

4 


0.7994912 


0.7992237 


3 X 10"^ 


0.8568976 


0.8568895 


8 X 10-6 


1 
5 


0.8887595 


0.8887601 


6 X 10"^ 


1 
6 


0.9090155 


0.9090167 


1 X 10"*^ 



Table 1: Numerical data compared with the analytic prediction for different values of the filling 
fraction a of operator Oi . 

which agrees neatly with the strong couphng result (32) within the numerical error. In figure 
6 we plotted several other fits and data for several different filling fractions. We see that, very 
non-trivially, the weak coupling results approach all the predictions from strong coupling in 
the Frolov-Tsetytlin limit! 

In appendix D we provide another simple example, still using the folded string as our 
toy model but using a different BPS operator = Tt{Z^) which has been recently studied 
abundantly in the literature [5, 6, 22, 7, 38, 39]. 

A few words of comfort to the readers less familiar with numerical fits. Throughout the 
paper we will perform several fits in order to extrapolate our finite J data and find the large 
J asymptotics. We are given a set of data {J, r(J)} with n points and perform a fit of this 
data with n — 1 coefficients, e.g. something like r = ao + ai/ J + • • • + an-i/ J^~^- We use the 
data to fix the constants ak using Mathematica. Then we plot the resulting fit for values of 
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J which are typicaUy much larger that the J's in the data (that is the whole idea of using 
fits!). The point to emphasize is that the only assumption we need to make is on the form 
of the expansion of r(J). Otherwise there is no fine-tuning whatsoever involved. 

3.2 Landau-Lifshitz coherent states 

In this section we shall study the classical limit of SU{?>) spin chains, see e.g. [21]. At 
each site we have three degrees of freedom which we denote as \X) , |y) and \ Z) . We use the 
notation |u) = Ui\X) -\-U2\Y) +Us\Z) to denote a linear superposition of these three degrees 
of freedom. We normalize u to be a unitary vector so that automatically 

(u|u) = u-u= 1. (34) 

The Hamiltonian we consider is the one arising in the one-loop computation of the Dilatation 
operator in A/" = 4 SYM. It reads 

A ^ 

^ = ^E(Wi-Pn,n+i) (35) 

n=l 

where Pn,n+i is the permutation operator which swaps the degrees of freedom in sites n and 
n + 1 while In,n+i is the identity. Generic eigenstates are Bethe eigenstates which generalize 
(3). They are very entangled complex quantum states. A huge simplification occurs in the 
large J classical limit mentioned in the introduction. In this limit the exact states can be 
approximated by coherent classical states which simply read [19] 

|(^) = |u(^))®...®|u(^)). (36) 

Note that the norm of such states is automatically equal to 1. In the operator language this 
state translates into (6). 

The physical picture is that of a slow varying spin wave pointing in direction u which 
slowly changes as we go from one site to the next. It is then convenient to package all the 
information about the classical state in a continuous field 

u(a) = u(^) with a = 271 J. (37) 

To get some experience with these classical states let us review the computation of the energy 

^^For example, the a = 1/3 curve in figure 6 has quite a nice non-trivial behavior with a minimum and an 
inflection point nicely converging from below to the predicted analytical result. This is the honest outcome 
of a fit like the one described in the text, no fine-tuning is involved. Similarly, in appendix D we see that 
the fits predict nice constant asymptotics for r(J) which are approached well beyond the range of the data, 
figure 10. The match with the analytics is again perfect without any fine-tuning whatsoever. 
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of such low wave length states 



A 

8^ 



J 



{^\H\^) 



^{nn\® (u„+l| (In,n+1 " Pn.n+l) |Un) ® |u„+l) 



n=l 



A 

8^ 



J 



u„ • u„+i)(u„+i • u„)] 



(38) 



A 

2J 




u . 



This expression precisely agrees with the strong couphng result (12). By identifying H\lp) = 
idr\^) in the left hand side of (38), we reproduce the integrated strong coupling equation 
of motion (19). More generally, the coherent states evolve in time as follows from the 
Schrodinger equation. That evolution is equivalent to the equation of motion (19) [19]. 

The Hamiltonian is just one of the many conserved charges of the integrable spin chain, 
which is a classical operator in the classical limit. The (trace of the) transfer matrix encodes 
all of them. Hence, a nice way of identifying whether a coherent state mimics an exact state 
is by measuring the expectation value of the transfer matrix and imposing that it matches 
the classical limit of the quantum transfer matrix [12]. 

A related but more technical comment is the following. The exact Bethe eigenstates 
are highest weights, S^\il)) = 0. The coherent states which we use to mimic them are only 
highest weights at the level of averages, that is {lp\S^\lp) = 0. Similarly, the exact states are 
cyclic, e'^^lip) = {ip) while {Lp\e'^^\Lp) = 1. Cyclicity of the exact states means that we don't 
need to care about where we break them and so on when computing expectation values or 
three-point functions. On the other hand, to mimic these computations with coherent states 
we should align the operators as done in the following section. 

3.3 Heavy-Heavy-Light three-point functions 

In this section we present the idea of the derivation of the general match between weak 
and strong coupling using the coherent state machinery reviewed in the previous section. 
A rigorous proof would involve cluttering the text with lots of kets, bras, tensor products 
etc without great benefit of clarity. Instead we sketch the derivation and the main physical 
ideas, leaving to the diligent reader the pleasure of filling in the details. 

One of the main issues in a rigorous derivation is what we mean by O2 been roughly the 
complex conjugate of Oi. Naively one might think that the precise details of how O2 differs 
from Oi should not matter. However we will show that it does and hence the complication, 
see sections 3.5 and 3.6. In any case, since we will be more precise about this point later, we 
will sacrifice rigor in this section. 

Having the necessary machinery, we move to the main object of study of this paper. 
We will consider the computation of three-point functions where two of the operators Oi 
and O2 are classical operators (well approximated by coherent operators/states) while the 
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third operator Os is a small operator. To simplify the expressions we will use the notation 
Ci = (length of operator Oi)/27r. 

For now, we take Os to be a vacuum descendant, 

Os=Af (Tr [X^^Y^'Z^'] + permutations) , 

where the normalization coefficient Af is simply one over the square root of the number 
of permutations. We only distinguish permutations up to a cyclic permutation due to the 
cyclicity of the trace 

To accomplish our businesses with combinatorics let us notice that among all the terms in 
Os only the ones of the form Tt[Z^^ anything] can give a nonzero contribution to the Wick 
contractions with Oi and (92, see figure 2b. The number of such terms is obviously equal to 

The operators Oi and O2 are described by coherent states as introduced in the previous 
section. In the operator language they are given by (6), while in the spin language we have 



= u(^)^® u(^gi)^®... (40) 

= •••®(u(^) 0(u{^) (41) 



The three-point function is now obtained by Wick contracting the three operators as depicted 
in figure 2b. Since there is no entanglement in operators Oi and O2 the contractions are 
trivial. The operator Os is glued, that is Wick contracted, with the other states at sites 
/c, /c — 1, ... . We should then sum over the insertion starting point k. The Wick contractions 
between Oi and O2 give 

L+k 

Ik= H u(i/£2) • u(i/A) . (42) 

i=k+l 

where L is the number of contractions between Oi and What is left is to consider the 
Wick contractions between and 02- Consider one of the terms in which gives 
a nonzero contribution Tt[Z^^XXY . . .]. The Wick contraction of this operator with the 
coherent states Oi, O2 at position /c, — 1, . . . gives 

^3(4)«3(y ) . . . ■ ■ ■ ■ 

Since the field changes slowly over those sites, we can approximate this contribution by 

Thus, all nonzero Wick contractions between and O2 give approximately the same 
contribution and there are {^^^^^^) of them. Hence, the result of these contractions is 
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All together, we get 

Ci23^J]X,J,, (43) 

Note that we have a gauge invariance when representing the coherent states as 

tensor products of single site states. That is, we can multiply the state associated with any 
site of the chain by an independent phase factor 

u07£) ^ e^^(^) u(i/£) , uOy£) ^ e-^(^) u(i/£) . (44) 

All observables, including the structure constant (43) are of course invariant under the gauge 
transformation however X]^ and J]^ in (43) do transform nontrivially (for ^ ^1+^2)- We can 
profit from this fact to choose the most convenient gauge. Namely we can fix a conformal-like 
gauge by demanding that Xk does not depend on k. That is 

^ u(fc//:2 + L/A) • u(fc//:i + L/A) 

Xk-i u{k/C2) • ^{k/Ci) 

We have 

nik/Ci + L/Ci) ~ u(fc/A) - ^-^^dMk/Ci) + Oil/Cj) 

and similarly for u. Thus we have to require that 

{ji+j2-js) u.a,u = 0, (45) 

which matches nicely the Virasoro constraint (20). In the conformal-like gauge (45), we 
further have that 

Jo = n u(^/£2) ■ ui^/C^) ~ exp J log (u • 1 - ^i±|^aa, u) = 1 + 0(1/ C) 

i=l Q 

We conclude that in this gauge: 

I '. 

^••-^V^-^Z-'S''-''"'"' '''' 

^ 

which is precisely the strong coupling result (25)! This is our main result. Notice also that 
if we substitute Va from eq.(9) instead of Ua we get 



which is indeed the hmit J a S> ja of the known BPS formula eq.(62). 
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3.4 Non-BPS operator O3 



Our tree-level computation did not rely in any crucial way on the small operator being 
protected. We can therefore trivially generalize it to any small SU{3) operator as 
we will now explain. All the manipulations of the previous subsection involving the large 
operators were not sensitive to the precise form of the small operator and only depended 
on its global charges. The only difference when considering non-protected small operators 
is in the combinatorial factor outside the integral in (46). That is, the more general weak 
coupling result is 

=^ j < ' (48) 



where T) is the sum of the coefficients of the traces in O3 in which all the Z fields are next 
to each other. 

To illustrate this point let us consider a simple example where is the Konishi operator 



1 



1 C-i 90 1/3 



T,{Z'X')-T,{ZXZX) ^ V=^ ^ 7^ = ^\h- (49) 



V3I ' ^ ' n V3 V3V2 



This ratio is exact and does not rely on the classical limit. Note that although in this 
example Os is in the SU{2) sub-sector the operators Oi and O2 can be generic SU{?>) 
operators. Another simple example is 



Tt{ZXY) -Tv{ZYX) V = 0. (50) 



If we restrict to the SU{2) sector considered in [4], see figure 1, we have^^ 



(51) 



where we have used that the dependence on Oi and O2 has canceled out and ^ were 
introduced in [4]. Recall that we are now considering non-BPS operators which can be 
parametrized by its Bethe roots; those are the Wj. Again, this formula is exact and does not 
rely on any classical limit. 

For non-protected operators however, we don't expect a match with strong coupling in the 
classical limit considered so far. That can be seen already at one loop. A main ingredient 
for the match between strong and weak coupling in the Frolov-Tseytlin limit is that the 
perturbative expansion organized itself in powers of X/J^ = I/k^. At one loop, there are 
two types of corrections [1]. One is by insertion of the one-loop Hamiltonian between the 
operators (see figure 9). The other is the two-loop correction to the Bethe wave functions. In 
contradistinction to the case where the small operator is BPS, these two types of corrections 
are not suppressed by and therefore, do not follow the Frolov-Tseytlin scaling. 



^^Three-point functions of two heavy operators and certain hght non-BPS operators were discussed in [7]. 
^^For in the 5/7(3) case, the first equal sign in (51) still holds provided we take the classical hmit. 
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If, on the other hand, we take the double scahng hmit 1 <C j <C J with A/j^ smaU then 
the operator O3 can also be treated classicaUy while its back-reaction on the large operators 
can still be neglected. In such limit, the weak and strong coupling structure constants 
might match. In the classical limit, we have computed the ratio ^ (51) in the past by 
simplifying the exact expression and reported the result in [4]. Combined with (51), we get^^ 



C--; = exp {jV + 0{f)) 



(52) 



where 



r = -iog 
J 



27r 



da 
27r 



1 ^3 



dt 



du 

2^1 



glog(2sin(tg/2)) — dup log (2 sinh(7rtp)) 



UCk 



Let us summarize the meaning of the several symbols in this formula. The operator O2 is a 
classical operator that can be approximated by an SU (2) coherent state parametrized by the 
Landau-Lifshitz field Ui{a) and Us{(j) as reviewed in the previous sections. The operator Oi 
is roughly the complex conjugate of O2 in the sense (C) discussed in section 3. The operator 
Os contains ji scalars X, js scalars Z and j = ji + js- The ji ^ 1 Bethe roots of this 
operator organize themselves into several cuts Ck- The density of roots in those cuts is p(v). 
Finally, 
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dv 



P(v) 



(53) 



UCk 



is a sort of trimmed quasi-momenta. For more details, we refer the reader to [12, 4, 40]. It 
would be very nice if that result could be reproduced at strong coupling. A first place to 
look at might be the small spike solution considered in [5]. 



3.5 Back-reaction numerically 

In this section we begin our study of back-reaction. As before, we start by some numerics to 
gain intuition about the general picture. We will again use the folded string introduced in 
section 3.1 as our toy model. A general analytical discussion is presented in the next section. 

As explained in the introduction the states or the eigenfunctions of the dilatation operator 
at one loop are labeled by the set of Bethe roots Uj = | cot y . Then the states corresponding 
to a classical finite-gap solution correspond to configurations of Bethe roots like the one shown 
in figure 4, where the roots are forming some dense cuts in the complex plane. 

The state O2 in general should have a diflFerent number of Bethe roots than the state Oi 
since the number of roots is related to the R-charge. Then, as discussed in the section 3.1, 
there are three major possibilities as to where the extra roots can go: 

(A) One option is to put extra roots at infinity (i.e. with zero momenta pj). Roots at infinity 
correspond to acting with a global symmetry generator on the classical configuration, 
see e.g. [34, 26]. 

^^Up to a pre-exponent. That is, F = 0(1) and this formula does not capture 1/j finite-size corrections 
to r. In particular, the a integral in F can be replaced by a its saddle point value. 
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Figure 7: Numerical data for the folded string example for the case (A) against the analyt- 
ical prediction (dashed gray lines) for different values of the filling fraction of operator Oi: 
a = 1/3, 1/4, 1/5, 1/6^ from bottom to top. 

(B) Another option is to add roots at finite values outside the classical cuts. 

(C) The last option is to add more roots to the classical cuts which exist already. 

In section 3.1 the last possibility was examined numerically and a perfect match with the 
analytical result (46) was found. In this section we discuss the numerical results concerning 
the first two options. 

3.5.1 Numerical results for roots at infinity 

Here we study the option (A) from the list above. We consider Os to be a BPS operator 
with j /2 scalars X and j /2 scalars Z for even j. In this case the final operator Oi and O2 
have the same length, but Oi has n = j/2 more Bethe roots than the operator We 
consider the situation when the extra n roots are sent to infinity, which means that Oi is a 
descendant of 

For the simplest case n = 1 we can compute the result exactly. In appendix C.l we show 
that 

^ J_ / l-2a + 2/J 
0^23 ^^a(l-a + l/J) ■ ^ ^ 

This suggests that in general one should have 

^^23 - ^ (55) 



1^ 



^000 Jn/2 ' 

^123 ^ ^ 
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Indeed the numerics for the case n = 2 and a = 1/5 give 
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0.27036118 
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0.092991262 


0.068076848 


0.053258866 



Fitting the data by X^^^o % obtain a = —1.01, which is indeed consistent with eq.(55).^^ 
Notice that for exactly the same O2 and but in the situation when aU roots belong to 
the big classical cuts of we obtained a finite value r = 0.88876 by fitting the numerical 
data to infinite length (see table 1). Thus we see that back-reaction is very important in this 
case and the result is completely different compared to the option (C) studied in section 3.1. 

In the next section we give a simple reason for that behavior and draw the general criteria 
for back-reaction. 

3.5.2 Numerical results for finite roots 

To keep the zero momentum condition in Oi, the two extra finite roots that we add should 
have opposite mode numbers k and —k. The case A; = 1 is the one considered in section 3.1; 
it amounts to adding more roots to the already open classical cuts. The cases k = 2,3,... 
corresponds to the case (B) from the list above. 

The main question is whether we get the same result when adding the two extra roots 
with mode number k 1 compared to the case k = 1 considered in section 3.1 or not. 
To answer this question we use again the analytic results [4] and evaluate these expressions 
numerically. The results are shown in figure 8. We conclude that again the answer is no and 
that the result diflFers considerably. In other words, back-reaction is relevant. 

3.6 General criteria for back-reaction 

In this section we will give a qualitative explanation for the dependence of the result on the 
type of modification of the state \Oi) with respect to the state \02) which is necessary for 
R-charge conservation. As we discussed above there are three options (see figure 5). The 
case (A) we mainly discuss in appendix C. One can think about the case (A) as a limit of (B) 
when the mode number k of the extra roots goes to zero. In this section we focus therefore 
on cases (B) versus case (C). 

To understand better the diflFerence between these two cases, we first consider the example 
of a free massless scalar field (/)(cr, r) on the cylinder. A general classical solution is of the 
form:^^ 

0(cr, t) = ^An COS (n(cr + r) + 0^) . 

n 

The analog of the finite gap solution would be the situation when the sum has only finitely 
many terms. The amplitude An is related to the number of excitations inside a cut through a 

^^We also considered cases with more Bethe roots at infinity and found again evidence for (55) although 
the precision decreases substantially as we increase the number n of Bethe roots sent to infinity. 
^^For simplicity, we only consider the left-moving solutions. 
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Figure 8: Numerical data for the three-point function with operator Oi deformed by adding two 
finite roots with mode numbers k = ±3 (option (B)) divided by the same three-point function with 
operator Oi deformed by adding two roots to the already existing cuts with mode numbers k — ±.1 
(option (C)) for different filling fractions of operator Oi: a = 1/3, 1/4, 1/5, 1/6 (from bottom to 
top). We see that the ratio does not go to one. In other words, back-reaction is relevant. 

Bohr-Sommerfeld quantization condition, which in the present case coincides of course with 
the osciUator nA^ = cN^ where c ^ /i is a smaU constant and A^^ is the number of excitations 
at mode number n caUed filling number. We see that 



Suppose we start with an unperturbed solution consisting of a single large cut 

(/)o(cr, 0) = Acos(na) . 

There are two different types of perturbations we can introduce. One is by adding a small 
number of excitations, /, to the existing large cut 



The other, is by adding the small number of excitations / at a different mode number k ^ n 

0) = Acos(na) + 5A^^^ cos(fca) . 

These perturbations obviously correspond to the options (C) and (B) in the interacting case. 
Since the amplitude is proportional to the square root of the filling number (56), for the case 
> / we get 




(56) 



cos(ncr) . 





(57) 
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We see that in these cases the scahng of the correction to the wave function is very different. 

Having understood this let us turn back to our calculation of the three-point function at 
weak coupling and see how the difference between the two behaviors in (57) really matters. 
The only difference between the quasi-classical approximation for the Bethe eigenstates and 
the free field example is that now the specific form of the solutions are more complicated and 
when we add new roots, the old roots get a small correction. These differences are irrelevant 
for the behavior (57). 

In the large J limit, the quasi-classical approximation is valid and we can repeat our 
previous calculation of the three-point function as prescribed in section 3.3. For simplicity 
we consider the case ^3 = ji + ^2 so that the lengths Ci = JC2 are equal and we do not have 
to worry about the U{1) gauge transformations (44). This time however, let us be a bit 
more careful about back-reaction. Namely we assume that the coherent states for Oi and 
O2 are a bit different and are parametrized by two different functions U2(cr) = u(a) and 
Ui(a) = u(cr) + 6u{a). Since Ui is a unit vector, we have 

ti • 5u + 5u • u = —6u • 6u . (58) 

An important step in the previous 'naive' calculation was the proof that the direct Wick 
contractions between Oi and O2 (42) are trivial. The overlap of the two coherent wave 
functions is equal to 



I = exp ( ^ / da log{l + u • 5u) \ 0^ exp (-^jda 



u - 6u- ^(u • Suf 



(59) 



As we are only interested in the absolute of the three-point function, we only need to consider 
the real part of the integrand. The real part of the first term in the integrand is precisely the 
left-hand side of (58). Thus we see that the real part of the integrand is of order ^ J||5u|p. 

For the type of correction (C) where the roots are added to the big open cuts, we found 
that 5u(^) - //J. In that case, the integrand goes as /^/J and the overlap I is indeed 1. 
That lead us to the main result (46). For the second type of corrections where the roots are 
not added to the big cuts, we found that ^u^^^ ^ ^//J. In that case we then get 

\I\ ^ e-^^ (60) 

and thus for this case we have an additional exponential suppression for large / (the coefficient 
must be negative since the overlap cannot exceed 1). Furthermore, in this case, (60) is not 
even the only correction. This is indeed the case according to what we found from numerics. 
That is, for small /'s we found an order 1 correction compared to case (C) studied in section 
3.1. 

Note that if we try to pass to the limit /c ^ corresponding to the case (A), the wave 
function correction goes to infinity as one can easily see from eq.(57). This indicates that 
when adding roots at infinity, the suppression should be even stronger in agreement with our 
finding (55). 

We end this section by noting that at strong coupling, the classical operators are described 
by the same coherent states. Our findings in this paper indicate that also at strong coupling, 
the three-point function of two large and one small operators is dominated by the overlap of 
the wave functions. As such, we expect the above criteria for back-reaction to apply. 
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4 Conclusion and open problems 



In this paper we considered the computation of structure constants (7i23 of single trace 
gauge-invariant operators in A/" = 4 SYM. 

We further focused on a classical limit where the operators Oi and O2 are made out of a 
large number of constituent fundamental fields while Os is a protected BPS operator which 
is taken to be small compared to the other two. This means that O2 is a classical operator, 
Os is a quantum operator and Oi is roughly the complex conjugate of 02-^^ Under a proper 
identification of Oi 0^ O2 we obtained a perfect match between the leading weak coupling 
results and the strong coupling results in the Frolov-Tseytlin limit, see also figure 3. The 
match holds for any classical string solutions with nontrivial motion in and which admit 
such limit. Our result is analogous to the match observed in the spectrum problem [12, 24] 
and we hope it can be equally inspiring in guessing an all-loop interpolation for structure 
constants as done for the spectrum in [25, 26]. 

At the same time we also found that the precise form of 0\ 0^ O2 does matter. We 
get a match for a particular (albeit natural) choice of Oi. This choice can be described as 
follows. The classical operator O2 is described by many Bet he roots which condense into 
cuts representing a classical wave function [12]. One such example is depicted in figure 4. 
Then, in the setup considered in this paper, Oi corresponds to an operator with some more 
Bethe roots. The question is then, where do we put them? The answer is quite simple: we 
should add them to the classical cuts that are already present in the classical operator 02-^^ 

The three-point function of two classical operators and one quantum operator is domi- 
nated by the overlap of the two classical operators. If instead of adding the extra roots to 
the open cuts we put them outside, then the change in the classical operator is not small 
enough. Putting roots outside the cuts at finite values can be interpreted as considering 
quantum fiuctuations around the classical solution [35] . Putting extra roots at infinity corre- 
sponds to acting with global symmetry generators on the classical solution. In both of these 
cases, the classical state Oi still seems to be very similar to O2. However, for these cases 
the structure constants are very sensitive to the precise details of the deformations and do 
not match the strong coupling result. The match only occurs when the roots are added to 
previously open classical cuts. 

In addition we computed the structure constant between three non-protected classical 
operators with one of the operators still much smaller than the other two, see (52). It would 
be great if this could be matched with strong coupling. The asymptotic analysis of the small 
spike from [5] might be of great use in this case. 

One of the main ingredients of this paper was the use of coherent states for studying 
the weak coupling classical operators. The other main ingredients were the exact tree-level 
formulae for the structure constants derived in [4] and evaluated numerically in the present 
paper. They allowed us to be confident about the validity of the coherent state approach 
and they also provided very convincing evidence for the physical picture vindicated in the 

^^We cannot have simply d = O2 because of R-charge conservation. 
According to our discussion in section 3.6, we expect that if the classical solution contains many classical 
cuts we can add the extra roots to any of the existing cuts and get the same classical result. 
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previous paragraphs. 

For example, for the folded string we solved the Bet he equations numerically and plugged 
the solution into our analytic results [4] . We then compared the answer with the analytical 
predictions from coherent states and found a precise agreement in the classical limit. 

Note that the exact expression for the type of structure constants considered here is very 
complicated [4]. Here, in the classical limit, we have obtained a huge simplification of it, see 
(46). It would be very interesting to understand how to derive the simplified result directly 
from the exact expression. Such understanding may teach us how to do the same for the 
case where all three operators are generical classical operators. 

For the computation of the classical three-point functions, all we really needed was the 
fact that the large operators are well approximated by coherent states. Such semiclassical 
coherent state approximation exists for any spin chain in the appropriate limit and hence 
our methods can be used for any large gauge theory. No need for integrability (the match 
with strong coupling is of course very special to A/" = 4 SYM). 

There are several interesting possible extensions of this work. One of these is to extend 
the analysis done in this paper for the SU{3) sector into other sectors as well as the full 
theory. At weak coupling, even for the full SU{3) sector the exact expression is not yet 
known. However, as done in this paper, it might be easier to derive the classical results 
directly in the coherent states language rather than to simplify the exact expressions. 

Another important direction would be to extend this work to one loop. It is simple to 
expand the strong coupling result (22) to next order in At weak coupling, one should 

decorate the tree-level contractions by inserting the one-loop dilatation operator [1, 2, 3] as 
represented in figure 9 and correct the wave functions. An optimistic scenario would be an 
agreement up to two loops. 

Finally, it would be interesting to extend our weak coupling results to higher-point func- 
tions. In some specific limits, we expect the same techniques to apply. One such case is 
the four-point function of four large operators that are designed such that their four-point 
function is dominated by the exchange of a finite set of small operators. Another case is the 
four-point function of two large classical operators and two small BPS operators [7, 39]. In 
this case we expect that only a simple modification of our result is needed. 

Of course, the most exciting future direction would be to address the question of the 
discretization of the classical results and then conjecture the all-loop result for the structure 
constants, at least for operators that are asymptotically large. In the past, for the spectrum 
problem, this line of attack was extremelly fruitful, see e.g. [25, 26, 23]. 

^^Note that the exact expressions for the two cases are written in terms of the same building blocks. In 
a similar case where only one of the three operators is classical whereas the other two are BPS, we have 
already succeed in simplifying the exact expression. The result in that case is an exponent of a non-trivial 
functional of the classical data [4]. 
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Figure 9: At one loop, the three-point function gets two types of corrections. One type is the 
two-loop correction to the external wave functions eigenstates. The other correction shown in this 
figure comes from the one-loop Hamiltonian insertion H between legs that originate from one of the 
operators and goes to two different operators. As opposed to the overlap of the wave functions, that 
contribution is very local and therefore it is somehow simpler. It can be thought of as measuring 
the energy cost of splitting the string [1]. Outside the scalar sector the picture seems to be more 
involved [3]. 
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A Tailoring numerics 

According to the results of [4] , the ratio C^^s / C'm can be computed in Mathematica as 

_ Cxxo[Ll,Nl,L2,N2,L3,N3,us,vs] 
^~ Cooo[Ll,Nl,L2,N2,L3,N3] ' ^ ^ 
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where Cxxo and Cooo are defined in appendix D of [4] and us and vs refer to tlie roots of 
operators Oi and (^2, respectively. For example, if we wanted to reproduce the first point 
for a = 1/4 in figure 6, we would simply run 

us={l. 98066463867 - 0.63872734235 1,1.98066463867 + 0.63872734235 I, 
-1.98066463867 + 0.63872734235 1,-1.98066463867 - 0.63872734235 I}; 
vs={2 . 35231505474 , -2 . 35231505474} ; 

r=Abs(§Cxxo [16 , 4 , 16 , 2 , 4 , 2 , us , vs] /Cooo [16,4,16,2,4,2] //FullForm 



to obtain r = 0.976868082298. 

In table 2 and appendix D of [4], the scalar product between part of Oi and part of 
O2 needed to compute Cxxo was written as a sum over all possible partitions of the roots 
involved (see equation (114) of [4]). However, from a computational point of view it is much 
more convenient to write this scalar product using the new recursion relation for SU (2) scalar 
products that we derived in equations (73), (77) and (78) of [4]. This is in fact what we did 
to obtain the numerical results of this paper. Despite this improvement, we were only able 
to consider cases where Oi had at most 14 roots. For example, it took about 10 hours on a 
quad-core processor to compute the last point for a = 1/3 in figure 6. 



B BPS 3-point function for SU{3) sector 

For the setup in (4) the structure constant between three BPS operators reads [13] 



Qooo \Ji+J2-Jl-j2/ V Jl-31 / V Jl / V ^ (62) 

' J \ (J-J?>\ ( j \ (h+j2\ ( J-jl-j2+j3 \ (Jl+J2-jl-j2\ 
^-JsA J2 Aji+JsA ji Aji+J2-ji-j2A Ji-ji J 



In the classical limit considered in this paper this expression simplifies to (24). 



C More on back-reaction 



In this appendix we will consider the case where a root is added at infinity and will show 
analytically that back-reaction is very important. We will start in section C.l with a case 
where (7*23 can be computed exactly without referring to coherent states. Then, in section 
C.2 we will study in more generality the infinite root problem. 

^^For simplicity, the discussion in this appendix refers to the case when all roots in the heavy operators 
are finite. After some slight modifications described in the main text of [4], the codes given in the appendices 
of [4] can (and were) used to compute the cases when some roots are at infinity. 
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C.l A special BPS operator O3 = Tr (ZX) 



In this section we will consider a specific example where (7*23 can be computed exactly. We 
will then use it to draw some general conclusions about back-reaction. The case in mind is 

Os = Tr {ZX) 

and for simplicity, we take Oi and O2 to be in the SU (2) sector. The point is that the Wick 
contraction of with Oi fiips one of the X fields into a Z field and therefore effectively 
acts as 5^+. Now suppose we take Oi to be a Bethe eigenstate. Since Bethe eigenstates are 
highest weight states, they are annihilated by 5^+, so C^^s = 0. 

Now suppose that, instead, we take 

where O2 is represented by an exact Bethe eigenstate \ip2)- Here S~ fiips one Z field into an 
X and is normalized such that it preserves the unit norm of \ip2)- We get that^^ 



C*23 — a/ -^3 — -^1 + 2 . 

For BPS operators, (62) reduces to 

In total we find that in the classical limit C^i's/C'm scales as 1/aA7, in agreement with the 
numerical experiments (55). We see again that the result is very sensitive to what exactly 
one means by O2 — Oi. In figure 7 we plotted the numerical results using the formulae from 
[4] against the analytic prediction just derived. A perfect match is obtained. 



C.2 Root at infinity, more general scenario and any coupling 

We consider now a more general scenario which includes the setup of the previous subsection 
as a particular case and which is furthermore valid at any coupling. It illustrates very clearly 
the importance of back-reaction. 

Consider three SU{3) operators O2 and Os which are eigenstates of the exact dilata- 
tion operator. We consider a small operator Os and two large operator O2 and Oi O2. 

Suppose we now change the small BPS operator Os by adding to it one more X field and 
one more Z field. Schematically, 

ds = OsUXUZ . 

At the same time we keep the large operator O2 fixed. Due to R-charge conservation, the 
the i?-charge of Oi must be changed accordingly. We can do that in many diflFerent ways by 

'V^s - Ji + 2Cr2*3° = {HS^S-\^2) = (H [S^. S-] \^2) = {HSzm = {Js - Ji + 2) 
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adding one more root to Oi. If back-reaction would have been negligible, the result would 
not depend on where we add that extra root and the new three-point function would have 
been of the same order as the original one. As by now is clear from the examples in the 
main text, that is not true. To see that in this general consideration, we choose to add 
the new root at infinity. That is, Oi = Af[S]^^ Oi], where S]^ is the i?-charge rotation that 
rotates Z field into X scalars and A/" = 1/ V Js ~ Ji is a normalization factor of order 1 / V^, 
constructed such that the two point function of the new operator is still normalized to 1. 
Let us now consider the ratio 

We will argue that R is of order 1/VJ, which means that, as in the folded string example 
above, a small change in Oi and Os leads to a different large J scaling and back-reaction is 
very important. 

The three point function is invariant under a global R-charge rotation of all operators. 
Hence we can trade the rotation of Oi by an action on Os and The action on O2 gives 
zero since it is an highest weight. The action on yields a new small operator. On the 
other hand, the normalization factor JV = 0(1/VJ) remains and hence R = 0(1/VJ) as 
advocated. 

A similar reasoning with more roots at infinity would lead to (55). 

D Os = Tr(Z^) and (ignoring) mixing with double traces 

In the main text we focused on operators O3 with ji + ^2 7^ scalars X, F , and js 7^ scalars 
Z. The reason for considering this setup was pointed out in [4]: in this setup all operators 
are Wick contracted to each other as depicted in figures 1 and 2. When this is the case, a 
simple large counting shows us that we can ignore mixing of the single trace operators 
with double traces. 

Another setup corresponds to ji = ^2 = 0, 

Os = Tt{Z') . (63) 

In this case, the length of operator O2 is the sum of the lengths of the other two operators. 
In other words, there are now no tree-level Wick contractions between Oi and Os. In this 
case to compute the structure constants we must also consider the mixing of operator O2 
with double traces [14].^^ The setup in [4] was chosen to avoid this complication. 

On the other hand, the case = Tt{Z^) is the case which has been most extensively 
studied in the literature recently at strong coupling [5, 6, 22, 7, 38, 39]. In all these computa- 
tions, the mixing with multi-string states has not been taken into account. Let us ignore the 
mixing with double trace effect and see what we get in the classical limit at weak coupling. 
Basically, the hope is that the same effect is being forgotten at weak and strong coupling 

^^At the same time, we should mention that there is however considerable evidence in the literature that 
extremal correlators can be considered as an analytic continuation of non-extremal ones [41], see also [42]. 
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Figure 10: Fits of the numerical data (black curves) compared to the analytical prediction (dashed 
gray lines) for different values of the filling fraction of operator Oi: a = 1/3, 1/4, 1/5, 1/6 for j = 3. 
No fine tuning is involved, see discussion at the end of section 3.1. 

and hence "cancels out" leaving a remainder quantity behind which can still be matched. 
Another plausible option is that this effect is negligible in the classical large L limit 

This seems indeed to be the case. From the point of view of the coherent state derivation, 
presented in section 3.3, we can freely set ji = j2 = 0. Therefore, the weak coupling result 
(46) matches again the strong coupling result in the Frolov-Tseytlin limit, (25). 

Now, let us be more precise and face the kind of question which we already encountered 
in other examples described in the main text, see e.g. section 3.1. Namely, what exactly 
are we computing? Given a classical operator (^2, what is the operator Oi Operator 
Oi has Ji scalars X and total length J. Operator O2 has Ji scalars X but a slightly larger 
length equal to J + j. Operator O2 is parametrized by some configuration of Bethe roots like 
the one in figure 4. The position of the roots is uniquely fixed by a choice of mode numbers 
and by the length of the operator. Hence, a most natural choice for Oi is to take it to be the 
(conjugate of) operator whose Bethe roots have the same mode numbers but slightly smaller 
length. 

To be more concrete, consider O2 to be the operator dual to the folded string introduced 
in section 3.1. Then Oi would be the same folded string with slightly smaller length. From 
the string point of view this would be a folded string with the same angular momentum in 
the plane identified by the scalar X and with j less units of angular momentum in the plane 
identified by the scalar Z. 

We performed once more a numerical check of the coherent state prediction (46) =(25) 
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against the results of [4] when we put ji = in the formulae of that paper. We find again a 
perfect match between the fits of the numerical data and the analytic prediction, as depicted 
in figure 10. 
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